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A Distributed and Privacy-Aware Speed Advisory 
System for Optimising Conventional and Electric 

Vehicles Networks 

Mingming Liu, Rodrigo H. Ordonez-Hurtado, Fabian Wirth, Yingqi Gu, Emanuele Crisostomi and Robert Shorten 


Abstract —One of the key ideas to make Intelligent Trans¬ 
portation Systems (ITS) work effectively is to deploy advanced 
communication and cooperative control technologies among the 
vehicles and road infrastructures. In this spirit, we propose a 
consensus-hased distributed speed advisory system that optimally 
determines a recommended common speed for a given area in 
order that the group emissions, or group battery consumptions, 
are minimised. Our algorithms achieve this in a privacy-aware 
manner; namely, individual vehicles do not reveal in-vehicle 
information to other vehicles or to Infrastructure. A mobility 
simulator is used to illustrate the efficacy of the algorithm, and 
hardware-in-the-loop tests involving a real vehicle are given to 
illnstrate nser acceptability and ease of the deployment. 

Index Terms —Speed advisory systems. Distributed algorithms. 
Optimisation. 


I. Introduction 

At present. Intelligent Speed Advisory (ISA) systems, as 
a part of Advanced Driver Assistance Systems (ADASs), 
have become a fundamental part of Intelligent Transportation 
Systems (ITS). Such systems offer many potential benefits, 
including improved vehicle and pedestrian safety, better 
utilisation of the road network, and reduced emissions. 
Recently, many papers have appeared on this topic 
reflecting the problem from the viewpoint of road operators, 
infrastmcture providers, and transportation solution providers 

Q-i)- 

In this paper, we consider the design of a speed advisory 
system (SAS) making use of vehicle-to-vehicle/infrastructure 
(V2X) technologies. Our starting point is the observation that 
different vehicle classes are designed to operate optimally at 
different vehicle speeds and at different loading conditions. 
Thus, a recommended speed, or speed limit, may be optimal 
for one vehicle and not for others. One possible way to 
handle this is by recommending a different speed to different 
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vehicles to take into account their differences (e.g., vehicle 
type, vehicle age, fuel mode, load, desired time of arrival). 
However, asking different cars to drive at different speeds 
does not make sense in any practical scenario. On the other 
hand, groups of cars following a common speed occurs 
frequently in practice. For example, roughly speaking, in 
highway driving cars tend to follow a given speed limit 
where possible. This tendency is increased in situations 
where intelligent tempomats are deployed. Other situations 
where common cars follow common speeds include dedicated 
lanes, and special zones in cities. Finally, generic benefits of 
following a common speed include, reduced emissions (due 
to less frequent accelerations/decelerations), reduced energy 
consumption, increased throughput, and increased safety 

0 - 113 - 

In this paper, we consider the problem of recommending a 
common speed to a set of cars. Our ISA will be constructed 
with an instantaneous optimality goal in mind. For example, 
we may wish to minimise the instantaneous group emissions, 
or the instantaneous group energy consumption. Specifically, 
we are interested in addressing the following questions: 

(i) Firstly, we are interested in designing an ISA system to 
recommend the same speed to a network of conventional 
vehicles (i.e., with an Internal Combustion Engine (ICE)) 
travelling along an extra-urban route (e.g., a highway), 
connected via a V2X communication system, such that 
the total emission is minimised; 

(ii) Secondly, we are interested in designing an ISA system 
to recommend the same speed to a network of EVs 
travelling in an urban network (e.g., the city centre), 
connected via a V2X communication system, such that 
the total battery consumption is minimised. 

We shall show that these problems can be solved in a 
privacy preserving manner without a large communication 
burden on the vehicle-to-vehicle (V2V) network. Eurthermore, 
we believe that algorithms addressing these problems will 
be of use when cars are asked to follow a constant speed 
to achieve certain goals. These include, transit through 
environmental and safety zones, highway driving with 
adaptive speed-limits (up to some maximum deviation from 
the nominal speed limit), and eco-driving for fleets of 
autonomous vehicles travelling in special lanes. 
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This paper is an extended version of the idea presented in 
m). In order to allow vehicles to collaboratively compute the 
optimal recommended speed, we shall further assume they are 
equipped with V2X technologies, can exchange information 
with their neighbours, and can exchange limited information 
with the infrastructure. We shall show that one can design, 
using very simple ideas, an effective SAS in a manner that 
preserves the privacy of individual vehicles. Extensive simu¬ 
lations, including hardware-in-the-loop (HIL) tests using real 
vehicles, are given to demonstrate the efficacy of our approach. 

II. Related work 

In this section, we give a brief review of some related 
work. First note that a detailed review of this topic is given 
in fT^ . Conventional systems are described in 0-0, HD, 
P6|. These papers describe various aspects of the ISA design 
process, which includes the design of driver display systems, 
the incorporation of external environmental information, and 
the algorithmic aspects of speed and distance recommenda¬ 
tions. Recently, there has been a strong trend to also include 
traffic density information; references describe 

work in this direction. In these works, density information is 
included in the procedure via loop detectors or via explicit 
density estimation using V2V technology. The differentiating 
feature of the approach followed in this paper is that density 
and composition of the vehicle fleet is also used, but in an 
implicit manner as part of the optimisation algorithm. Finally, 
we note that there is a huge body of work on cooperative 
control of vehicles and its connection to consensus algorithms 
|Tg, It is important to note that we are designing a 

SAS and not a cooperative control system. This distinction is 
important as it allows us to ignore string stability effects which 
are a fundamental limitation of many cooperative control 
architectures ||25l-||28l. 


recommendations are made. Let Si (k) be the recommended 
speed of the vehicle i G N := {l,2,...,iV} calculated at 
time instant k. Thus, the vector of recommended speeds for 
all vehicles is given by s (fc)^ := [si (fc), S 2 (k ),..., sjv (k)], 
where the superscript T represents the transposition of the 
vector. Note that between two consecutive time instants 
{k,k + 1), the recommended speeds are constant while the 
driving speeds are time-varying real-valued variables. We 
denote by Nl the set of neighbours of vehicle i at time instant 
k, i.e., those vehicles which can successfully broadcast their 
recommended speeds to vehicle i. 


In addition, we assume that each vehicle i can evaluate a 
function fi that determines its average emissions, in the ICE 
case study, or its average energy consumption in the EV case, 
were it to be travelling at the recommended speed Si {k). Such 
functions are typically convex functions of the vehicle speed, 
as it will be shown in more detail in the following sections. 
We shall further assume that these functions are continuously 
differentiable and with a Lipschitz continuous first derivative 
/' which is assumed with positive bounded growth rate, i.e.. 


0 < dZn < 


, ^max5 

a — 0 


( 1 ) 


for all a, 5 G M such that ay^b, and suitable positive constants 

(i) fi) 

'^min, timax. A Schematic diagram of the above is illustrated 
in Fig. [T] In this context, we consider the following problem. 


Problem 1: Design an ISA system for a network of vehicles 
connected via V2X communication systems, in order to 
recommend a common speed that minimises the amount of 
emissions, or the total energy consumption, of the whole fleet 
of vehicles. 


III. Model and Algorithm 
A. Problem Statement 

Our objective is to develop an algorithm that converges 
to a common recommended speed for a fleet of vehicles 
driving along the same stretch of a road (e.g., a highway), 
or in the same area (e.g., in the city centre). Let N denote 
the total number of vehicles in the fleet, and let us assume 
that they all are able to receive the ISA broadcast signal. 
For this purpose, each vehicle is equipped with a specific 
communication device (e.g., a mobile phone with access to 
WiFi/3G networks) so that they are able to receive/transmit 
messages from/to either nearby vehicles or available road 
infrastructure (e.g. a base station). We assume that each 
vehicle can communicate a limited amount of information 
with the infrastructure, the infrastructure can broadcast 
information to the entire network of cars, and each vehicle 
can send a broadcast signal to its neighbours. 

For convenience, we assume that all vehicles have 
access to a common clock (for example, a GPS clock). 
Let k G {1,2,3,...} be a discrete-time instant in which 
new information from vehicles is collected and new speed 


Base station 



Fig. 1; Schematic diagram of the framework |29 


The optimisation problem that needs to be solved in order 
to address Problem 1 can be formulated as follows; 


min } fj (s,) , 
jgN 

s.t. Si = Sj, fli j G N. 


( 2 ) 
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The above is an optimised consensus problem and can 
be solved in a variety of ways (for example using ADMM 


-|32|). Our focus in this present work is not to construct 
a fully distributed solution to this problem, but rather to 
construct a partially distributed solution which allows rapid 
convergence to the optimum, without requiring the vehicles 
to exchange information that reveals individual cost functions 
to other vehicles. This is the privacy preserving component 
of our problem statement. 


Comment; Note that (|^ involves the recommended speeds Sj 
and not the actual speeds at which the vehicles are travelling 
in reality. This implies, among other things, that the optimal 
solution of is the same, disregarding whether vehicles 
follow in reality to the recommended speed or not. 


Comment: In the following, we shall not solve © in a 
single step. Instead, we shall propose an iterative algorithm 
that at each step yields individual recommended speeds that 
will eventually converge to the same value on the consensus 
constraints. 


To solve (|^ we use the iterative feedback scheme 


s (fc + 1) = P (fc) s (fc) + G (s (k)) e, 


(3) 


pN 


where {P(fc)} € 
matrice^ e G I 
equal to 1, and G 
with some assumptions 


xN 


is a sequence of row-stochastic 
is a column vector with all entries 
I—R is a continuous function 
to satisfy 


as we 


shall 


see in 


Theorem Algorithms of this type were proposed and 
studied among others in p^-| 351; the principal theoretical 


contribution here is to extend this framework to a new class 
of optimisation problems and to give conditions guaranteeing 
their convergence. 


Remark: Eq. ([^ describes a consensus algorithm with an 
input. In the first part of the right-hand side of §, P{k) 
creates pressure on s{k) so that all components of s{k) 
achieve a common value. The second part of ([^ describes the 
regulation constraint that the nonlinear function of s(fc) must 
satisfy. Thus, we can use Eq. 0 to force all components 
of s{k) to achieve a common value and at the same time 
to satisfy a nonlinear constraint. The choice of G thus 
determines which constraint is satisfied when s{k) reaches 
the equilibrium. In principle, many choices of G can be used. 
In this paper, we choose G to solve an optimisation problem. 
The principal mathematical considerations in deploying this 
algorithm are to determine the mathematical conditions on 
G so that ([^ actually converges to the desired value. While 
this is not a primary concern of this paper, rather we are 
interested in the applications of this, in the remainder of this 
section we shall give a brief overview of such considerations. 


We will require that (|^ has a unique solution. Note that, it 
follows from elementary optimisation theory that if all the /^’s 
are strictly convex functions, then the optimisation problem (|^ 

’Square matrices with non-negative real entries, and rows summing to 1. 


has a solution if and only if there exists a y* G M satisfying 

N 

i=i 

In this case, y* is unique by strict convexity, and the unique 
optimal point of Q is given by 

s* := y*e G (5) 


In order to obtain convergence of ([^ we select a feedback 
signal 

N 

G(s(fc)) = -/iE/;(s,(fc)), (6) 

i=i 

and we obtain the dynamical system 


N 


s{k + l) = P{k)s{k)-y'^fj{sj{k))e, y G K. (7) 


i=i 


In 1^ it is shown that if {P{k)}^^^ is a uniformly strongly 
ergodic sequenc^ and p is chosen according to 

-1 


N 

V- ^U) 


( 8 ) 


0 < ^ < 2 , / 

\j=i 

then 0 is uniformly globally asymptotically stable at the 
unique optimal point s* = y*e of Q. Such systems were 
studied in and formally analysed in | [29| . Eor complete¬ 
ness, we formally state the relevant results from these works 
as a theorem (Theorem 1). An overview of the proof is given 


in the appendix, and the interested reader may refer to |29| 
for details. 


Theorem 1 ( | |29| |) Consider the optimisation problem 0 
and the optimisation algorithm 0, and the associated Lure 
system 

yik + 1) = h{y{k)), 

h{y) ■■= y + G{ye). 

If G is defined by 0 and the condition 0 holds, then the 
following assertions hold: 

(i) There exists a unique, globally asymptotically stable 
fixed point y* G R o/ the Lure system 0- 

(ii) The fixed point y* of (i) satisfies the optimality condi¬ 
tion 0. and thus y*e G R'^ is the unique optimal point 
for the optimisation problem 0- 

(Hi) If, in addition, {P{k)}j,^^ C R^^^ is a strongly 
ergodic sequence of row-stochastic matrices, then y*e is 
a globally asymptotically stable fixed point for system 0. 

An outline of the proof can be found in Appendix To 
apply the Theorem 1 to solve the optimisation problem we 
proceed as follows. Eor each k we define the P (k) as 

P^Ak)=\ ifj&Nl, , ( 10 ) 

[ 0, otherwise. 

^That is, for every /cq G N the sequence P(ko), P{kQ l)P (ko), .. 
P(ko +£)••• P{ko),. . . converges to a rank one matrix. See (29| for further 
details. 
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where i,j are the entries’ indexes of the matrix P (k), and 
r]j G M is a weighting factor. For example, a convenient choice 
G (^0, where 1*1 denotes cardinality, 

giving rise to an equal weight factor for all elements in the 
reference speed vector s (fc). 

The assumption of uniform strong ergodicity holds if the 
neighborhood graph associated with the problem has suitable 
connectedness properties. If sufficiently many cars travel on a 
given area, it is reasonable to expect that this graph is strongly 
connected at most time instances. Weaker assumptions are 
possible but we do not discuss them here for reasons of 
space; see for possible assumptions in this context. 

Now, we propose the Optimal Decentralised Consensus 
Algorithm for solving 0 as shown in Algorithm The 
underlying assumption here is that at all time instants all cars 
communicate their value /j {sj (fc)) to the base station, which 
reports the aggregate sum back to all cars. This is precisely 
the privacy preserving aspect of the algorithm, as cars do not 
have to reveal their cost functions to anyone. Also implicit 
information as derivatives of the cost function at certain speeds 
is only revealed to the base station but not to any other agent 
involved in the system. 


Algorithm 1 Optimal Decentralised Consensus Algorithm 

1 : for /c = 1,2,3,.. do 
2 : for each i G N do 

3: Get F {k) = X] f'j {sj (fc)) from the base station. 

jeN 

4: Get Sj (k) from all neighbours j G 

5: Do qi (k) =rii- Y, i^j (*) “ s* (k)). 

jeNi 

6: Do Si (fc + 1) = Si {k) + Qi (k) — fi ■ F {k). 

7: end for 

8 : end for 


Comment: We note that in any practical implementation of the 
previous algorithm, the recommended speed may be bounded 
above and below by the road operator (i.e., should the optimal 
solution computed by Algorithm exceed speed limits, or be 
unreasonably low). 

IV. A Speed Advisory System eor Conventional 
Vehicles 

In this section, we evaluate the performance of Algorithm 1 
using SUMO p7) . SUMO is a microscopic, open source road 
traffic simulator, which is frequently used for validation and 
prediction purposes in the ITS community. First, we evaluate 
it using conventional simulations, and then use a real vehicle 
embedded into a hardware-in-the-loop (HIT) emulation. In 
particular, we perform the following experimental activities: 

1) First, we compare the pollution emission of a fixed fleet 
of vehicles (with initial non-optimal speeds) before and 
after Algorithm 1 is turned on. 

2) We then provide simulation results which are dynamic 
in nature, where vehicles enter and leave a section of 


the highway where the SAS is active; in doing so, we 
provide a discussion about the rate of convergence of 
Algorithm 1, and also evaluate what happens when some 
vehicles decide not to follow the recommended speed. 

3) Finally, we give a HIL emulation with a real vehicle 
travelling on a real road, embedded into a emulated 
network with a fixed number of simulated vehicles. 

In all the above experiments, we evaluate the performance 
of Algorithm 1 for p = 0.001 and /r = 0.01 (a convergent 
set-up), and all the vehicles are also subject to the driving 
constraints arising from the SUMO simulation (e.g., in terms 
of acceleration/deceleration profiles). 


A. Cost funtions to represent pollution emission in ICE vehi¬ 
cles 

We shall adopt the average-speed model proposed in p8] to 
model each cost function fi in function of the average speed 
s as 

a-I-bsi-I-cs^-I-ds^-I-es^-I-fs®-I-gs® ^ 

where a, b, c, d, e, f, g, k G K are used to specify different 
levels of emissions by different classes of vehicles. In addition, 
we use the emission types from shown in Table |^and Fig. 
1^ corresponding to petrol cars/minibuses with up to 2.5 tons of 
gross vehicle mass. Also, we use the following vehicle types: 

• Type 1: accel. 2.15 m/s^, decel. 5.5 m/s^, length 4.54 m. 

• Type 2: accel. 1.22 m/s^, decel. 5.0 m/s^, length 4.51 m. 

• Type 3: accel. 1.75 m/s^, decel. 6.1 m/s^, length 4.45 m. 

• Type 4: accel. 2.45 m/s^, decel. 6.1 m/s^, length 4.48 m. 


TABLE I: Emission factors for some CO 2 emission types 
reported in |[^. Here, {e,f, g} = 0 and k = 1. 


Type 

a 

b 

c 

d 

R007 

2.2606E+3 

3.1583E+1 

2.9263E-1 

3.0199E-3 

ROM 

2.5324E+3 

6.8842E+1 

-4.3167E-I 

6.6776E-3 

R021 

3.7473E+3 

1.0571E+2 

-8.5270E-1 

1.0318E-2 

R040 

1.2988E+3 

2.0203E+2 

-1.5597E-0 

1.2264E-2 



Eig. 2: Curves for the CO 2 emission types in Table Black 
circles mark the minimum point on each curve. 
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B. SUMO Simulations with a Fixed Number of Vehicles 

In this experiment we consider 40 vehicles travelling along 
a highway. The set-up for this set of experiments is as follows. 

• Road; A straight 5 km long highway with 4 lanes. 

• Duration of the simulation; 600 s. 

• Algorithm sampling interval; AT = Is. 

• Switch-on time; the algorithm is activated at time 300 s. 

The experimental results are given in Fig. and Fig. As 
we can see from Fig. making a small change in a common 
non-optimal initial speed yields a significant improvement 
in CO 2 emissions, as about 172 g of CO 2 are saved every 
kilometer. If this result is integrated over a whole road 
network, and integrated over a longer time horizon (e.g., a 
year), signihcant improvements can be observed. For instance, 
the savings are comparable to those that can be achieved if 
a EURO-3 vehicle (e.g., R039) is substituted by a EURO-4 
vehicle (e.g., R040), where the savings are around 10 g/km. 

In the second experiment (see Eig. |^, we assume that the 
vehicles are initially driving at random speeds around the 
optimal one. Then, once the SAS is activated, the vehicles 
change their speed to eventually converge to the optimal 
one. Significant CO 2 savings can be noticed in this example 
again; however, carbon savings here are less than in the hrst 
experiment, since now some vehicles were already driving at 
(or very close to) the optimal speed. 


Evolution of the recommended speeds 



Total CO 2 emissions 



Eig. 3; Simulation results for the static case, before and after 
the activation of the algorithm at time step 300 s. Setup; 32 
vehicles of emission type R007 and 8 vehicles of emission 
type R021, with uniform distribution of vehicle types. 


C. SUMO Simulations with a Dynamic Number of Vehicles 

We now consider a dynamic scenario, with a time-varying 
number of vehicles. To do this we partition the highway into 


Evolution of the recommended speeds 


90 



200 250 300 350 400 450 500 550 600 

Time [s] 

Total CO2 emissions 



Eig. 4; Simulation results for the static case, before and after 
the activation of the algorithm at time step 300 s. Setup; 
uniform distribution of emission types in R014, R021 and 
R040, with uniform distribution of vehicle types. 

three consecutive sections LI, L2 and L3. We then proceed 
as follows. Eirst, vehicles enter the uncontrolled section LI, 
with constant speed (randomly chosen in a given range). 
After completing LI, vehicles enter the section L2 where the 
SAS is always active, so the vehicles iteratively update their 
recommended speed using Algorithm 1. After completing L2, 
they hnally enter section L3 where they travel freely. The 
experiments are setup as follows. 

• Road; three consecutive straight edges; 

- LI; 5 km long highway with 4 lanes, uncontrolled; 

- L2; 5 km long highway with 4 lanes, ISA controlled; 

- L3; 5 km long highway with 4 lanes, uncontrolled. 

• Total number of cars; 650, with uniform distribution of 
both emission types among R014, R021 and R040, and 
vehicle types. 

• Vehicular flow entering LI; one new car every 2 seconds 
until simulation time 1300 s. 

• Length of simulation; 3010 s. 

• Window size for the calculation of the moving average 
(MA) of CO 2 emissions for visualisation purposes; 500 
time steps. 

• Travelling speeds for cars on LI are randomly chosen 
with uniform distribution in 3 scenarios; 

- Case 1, constant speeds in (80,100) km/h. 

- Case 2, constant speeds in (60, 80) km/h. 

- Case 3, constant speeds in (40, 60) km/h. 

Note that even though this is a dynamic situation, the 
vehicle density on each part of the road becomes almost 
constant after certain time. A sample of simulation results 
is given in Eig. which reveals what might be expected 
from the initial experiments. Namely, the further vehicles are 
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away from the optimal speed, the more is to be gained by 
deploying the ISA. 


Moving average (MA) of CO2 emissions per section 



Fig. 5; Simulation results for the dynamic case: MA of CO 2 
emissions on LI and L2, for initial speeds on LI in Case 3. 


To generalise the previous results, we conducted one 
hundred random experiments for each of the three cases 
described above. In each experiment, we collected the 
simulation data of the total CO 2 emission generation on each 
section of the highway from SUMO. Table |n] summarises the 
aggregated results of this exercise, and clearly demonstrates 
the benefits of the ISA. 


TABLE II: General results: total emissions per lane. 



Total Emissions* [tonnes] 


LI (Uncontrolled) 

L2 (Controlled) 

Improvement 

Case 

Mean 

cr 

Mean 

u 

Mean 

cr 

1 

1.5199 

0.00114 

1.4684 

0.0002 

3.40% 

0.07 

2 

1.4751 

0.0004 

1.4648 

0.0001 

0.69% 

0.03 

3 

1.5945 

0.0028 

1.4679 

0.0001 

7.94% 

0.16 


*: Sum of emissions at every time step (i.e. time integration). Mean: average 
of 100 different measurements, cr: standard deviation. 

Comment; Note that the solution we have obtained is optimal 
for the environment and for the collective, e.g., in terms of the 
reduction of overall emissions. However, the solution might be 
unfair for some single users who would be recommended to 
drive at a different speed than originally desired. One way to 
improve fairness could be to decrease road taxes for virtuous 
vehicles, to compensate them from the inconvenience caused 
by the dirty vehicles in terms of recommended average speeds. 

Comment; Note that it is clearly the case that the design 
parameters of the algorithm have the potential to affect 
emission savings. For example, the speed of convergence of 
the algorithm affects the rate of which the emissions are saved. 

Following the previous comments, convergence issues are 
briefly discussed with the support of Fig. As we have 
mentioned in the previous paragraph, the recommended 
speed is updated in an iterative fashion, and it eventually 
converges to a steady-state value after a number of steps. The 
time of convergence is mainly affected by the topology of 
the communication graph and, in particular, by the number 


of neighbours. Convergence rate of algorithms of graphs 
is a well-studied topic. The interested readers are referred 
to | |3^ , pO] for the favours of results available. In this 
paper, we just give a brief empirical illustration. As we 
shall see, roughly speaking, the larger is the communication 
range of the vehicles, the larger will be the number of 
neighbouring vehicles, and (on average) the earlier is the 
convergence to the steady-state value. In particular. Fig. 
shows how convergence of a single vehicle is affected by its 
communication range. Note that convergence occurs after a 
few steps (less than three minutes) when the communication 
range is of the order of hundreds of meters. On the other side, 
convergence is very slow when the communication range is 
of only 15 meters, where it is frequent to get an empty set of 
neighbouring vehicles. Finally, observe that the steady-state 
value of the optimal speed is not constant as the optimal 
speed depends on the vehicles currently driving along the 
road section L2, and thus the types of vehicles (and in turn 
the optimal speed) change continuously. 


Evolution of the recommended speed 




Fig. 6: Evolution of the recommended speed (top) and number 
of neighbours (bottom) of a target vehicle, as functions of the 
radius r of its communication range. 


As has been previously remarked, convergence is for the 
recommended speeds, and thus it is independent on whether 
vehicles do follow the recommended speed or not. However, 
not following the recommended speed usually decreases the 
benefits of the proposed procedure in terms of CO 2 savings. 
Accordingly, Fig. illustrates how CO 2 savings are affected 
by the percentage of vehicles that decide not to follow the 
recommended speed. 

D. SUMO-based hardware-in-the-loop (HIL) emulation 

Finally, to give a feeling to a real driver of how this system 
might function, we now describe a HIL implementation of the 
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Total CO 2 emissions on L2 



Percent of vehicles ignoring the recommended speed [%] 

Fig. 7: Total CO 2 emissions as a function of the percentage 
of vehicles not following the recommendations; emissions are 
maximum when all vehicles do not follow the recommended 
speed, and minimum when all they do. Mean (x) and standard 
deviation (I) are calculated based on five different experiments. 


algorithm. Specifically, we use a SUMO-based HIL emulation 
platform that was developed at the Hamilton Institute ED’ 
p2| . This emulation platform uses SUMO to simulate a real 
environment and generate virtual cars, along with a dedicated 
communication architecture supported by TraCI (a Python 
script implementing a TPC-based client/server architecture) to 
provide on-line access to SUMO, a smartphone connected to 
the 3G network and running the plug-in SumoEmbed (designed 
for use with Torque Pro |43|), and a OBD-Il adaptor p) to 
embed a real car into the simulation, as shown in Fig. [^The 
idea then is to allow a person to drive a real vehicle on a real 
street circuit, to experience being connected to a network of 
(virtual) vehicles driving along a virtual environment based 
on the physical street circuit. Specifically, we performed this 
experiment by driving a Toyota Prius on a single-lane street 
circuit in the North Campus of the Maynooth University, while 
the Prius is embedded into the HIL emulation and represented 
by an avatar which interacts with the avatars of 29 other 
virtual (simulated) vehicles driving along the same stretch of 
(emulated) road. 



Real car 


Workstation 


Smartphone 

(SumoEmbed) 

■ 


OBD-II Adaptor 


Python 

script 


Server 

component 


Main component 


SUMO 


Fig. 8: Schematic of the SUMO-based HIL emulation plat¬ 
form. 


The experiment begins when the simulation is started on 
the workstation and the server component of the Python 
script waits for a call from the OBD-II connected smartphone 
in the real vehicle. Since the selected street circuit only has 


one lane, the vehicles are released sequentially from the same 
starting point. The avatar representing the Prius departs in 
the seventh position. Once the connection between the Prius 
and the workstation is established, the position and speed 
of the Prius’ avatar are updated using real-time information 
from the Prius via the OBD-II adaptor. From the point of 
view of the ISA algorithm, the Prius is regarded as a normal 
agent in the SUMO simulation, i.e. treated just like any other 
simulated vehicle. 

The consensus algorithm for the proposed ISA system is 
embedded in the main component Python script. Thus, once 
the respective recommended speeds are calculated, they are 
sent to the vehicles through the server component, via the 
cellular network to the smartphone in the case of the Prius, 
and via TraCI commands in the case of the other vehicles 
in the simulation. Note here that the driver behaviour is 
different for a simulated vehicle compared to the case of the 
Prius: while we force each simulated vehicle to follow the 
recommended speed as far as possiblE| the Prius’ driver is 
allowed to either follow or ignore the speed recommendation 
(displayed on the smartphone’s screen) as desired. 

The HIL experiment is setup as follows. 

• Length of the experiment; 600 s, of which the ISA 
algorithm is only engaged at around time 300 s; 

• Total number of cars: 30, with uniform distribution of 
emission types between R014 and R040, and uniform 
distribution of vehicle types, with a maximum speed of 
100 km/h. 

• The sampling time interval AT for collecting new 
information and updating the recommendations is 1 s. 

Results of the experiment are depicted in Fig|^ From Fig|^ 
it can be observed that the total emissions oscillate around an 
average value of 4323 g/km once all the 30 vehicles are added 
to the simulation, and also that they start reducing once the 
ISA algorithm is switched on, up to a hnal average value of 
4253 g/km. Finally, as can be observed in Fig|^ (bottom), the 
recommended speed can be roughly followed (on average) by 
the driver subject to the physical constraints on the real street 
circuit (e.g. traffic calming devices), and the discretisation step 
for the visualisation of the recommended speed. 

V. A Speed Advisory System for Electric Vehicles 

In this section, we extend our previous discussion to the 
special case of a fleet of EVs. The motivation for doing 
so is the recent interest in EVs as a cleaner alternative 
to conventional, more polluting, ICE vehicles. In recent 
years, some cities have decided to close the city centres to 
normal traffic, only allowing the transit to specihc categories 
of low (or zero) polluting vehicles; see for instance | [45) . 
Similarly, in many cities all around the world, the urban 
public transportation fleet has been restricted to electrically 

^Concerning mainly the interaction between vehicles and the design param- 
eters for the simulated cars such as acceleration, deceleration, car following 
model or driver information. 
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Fig. 9: Evolution of the overall CO 2 emissions (top) and of the 
speeds related with the Prius (bottom) for the HIL simulation. 
The algorithm was turned on around time 300 s. Average speed 
calculated with a window size of 20 time steps. 


driven vehicles, as for the case of fleets of urban electric 
buses; see for instance the cases of Sao Paulo in in Latin 
America, Louisville in the US ED’ or Wien in Europe pS) . 

Clearly, the previously designed Algorithm I needs to be 
adapted when applied to a fleet of EVs concerning different 
aspects. Eirst, in terms of the field of application, EVs are 
typically used for short distances due to their reduced driving 
ranges, and thus are most likely deployed in city centres; 
second, the cost functions should not consider instantaneous 
emissions that, in the case of EVs, can be considered equal to 
zero. Accordingly, in the following we show how the previous 
framework can be adapted to the case of a fleet of EVs, and 
now the objective is to maximise the energy efficiency of the 
fleet of cars or, in other words, to extend their driving range. 

A. Cost functions to represent energy consumption in EVs 

Most of the discussion here follows the reference p9l , 
where the ranges of EVs are reported for different brands and 
under different driving cycles. Power consumption in an EV 
driving at a steady-state speed (along a flat road) is caused 
mainly by four sources; 

• Aerodynamics power losses: they are proportional to 
the cube of the speed of the EV, and depend on other 
parameters typical of a single vehicle such as its frontal 
area and the drag coefficient (which in turn depends on 
the shape of the vehicle). 

• Drivetrain losses; they result from the process of con¬ 
verting energy in the battery into torque at the wheels 
of the car. Their computation is not simple, as losses 
might occur at different levels (in the inverter, in the 


induction motor, gears, etc); in some cases, these power 
losses have been modelled as a third-order polynomial, 
whose parameters have been obtained by fitting some 
experimental data (see EE)). 

• Tires: the power required to overcome the rolling distance 
depends on the weight of the vehicle (and thus, on the 
number of passengers as well), and is proportional to the 
speed of the vehicle. 

• Ancillary systems; this category includes all other 
electrical loads in the vehicle, such as HVAC systems, 
external lights, audio system, battery cooling systems, 
etc. Here, the power consumption does not depend on 
the speed of the vehicle and can be represented by a 
constant term that depends on external factors (e.g., 
weather conditions) and personal choices (desired indoor 
temperature, volume of the radio, etc). According to 
experimental evaluations | |49) , the power losses due 
to ancillary services usually vary between 0.2 and 2.2 kW. 


Thus, by summing up all the previous terms, the power 
consumption Peons can be represented as a function of the 
speed V as 


Ckg 2 

-h ai -I- a2V + asv , 

V 


( 12 ) 


where the left hand side is divided by the speed in order to 
obtain an indication of energy consumption per km, expressed 
in kWh/km. Such a unit of measurement is usually employed 
in energy-efficiency evaluations, and we shall assume that 
every single EV will use ( [T^ as its personal cost function f. 
Accordingly, Eig. 10a shows a possible relationship between 
speed and power consumption, obtained using data from 
Tesla Roadster and assuming a low power consumption for 
ancillary services of 0.56 kW (i.e., assuming air conditioning 
switched off). As can be noted from Eig. |10a| there is a large 
energy consumption at large speeds due to the fact that power 
increases with the cube of the speed for aerodynamic reasons; 
however, it is also large for low speeds, due to the fact that 
travel times increase and, accordingly, constant power required 
by ancillary services demands more energy than the same 
services delivered with high speeds. 


B. Experimental results 

According to the previous discussion, we now assume that 
the objective is to infer the optimal speed that the ISA system 
should broadcast to a fleet of EVs travelling in a given area 
of a city (e.g., in the city centre). Eor this purpose, we assume 
that a fleet of 100 vehicles travels in the city centre for an 
hour, and following the next steps: 

• In the first 20 minutes, the vehicles travel at the optimal 
speed calculated from Algorithm 1. 

• In the second 20 minutes, they travel at a speed below 
the optimal speed. 

• In the last 20 minutes, they travel at a speed above the 
optimal speed. 


In the first stage we assume that the communication graph 
among the EVs changes in a random way, i.e., at each time 
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Step an EV receives information from a subset of vehicles 
belonging to the fleet. This is a simplifying assumption that 
can be justified by assuming that in principle all vehicles 
might communicate to all the other vehicles (i.e., they are 
relatively close), but some communications might fail due 
to obstacles, shadowing effects, external noise, etc. Besides, 
in the two last stages we assume that the change of speed 
occurs almost instantaneously, since there is no requirement 
to iteratively compute an optimal speed. 


We tuned our parameters in Algorithm[^as p = /i = 0.001, 
and we simulate different cost functions for each EVs by 
assuming a random number of people inside each car (between 
1 and 5 people) with an average weight of 80 kg, and 
by assuming a different consumption from ancillary services 


while the average energy consumption is shown in Eig. 11b 


of the 


10b 

g- 

11a 


lb 




(b) All the cost functions overlapped. Black circles mark the minimum 
point on each curve. 

Eig. 10; Curves for the cost functions used in the experiment. 
All of them were chosen convex. 


As is shown in Eig. 1 lb a simple ISA can be used for the 
case of a fleet of EVs just as in the case of conventional cars 
(for what regards the mathematical background) by simply 
adapting the cost functions to the new application of interest. 


VI. Conclusion 

In this paper we present a new ISA system. The system 
is based on a solution to an optimised consensus problem. 
We show that the ISA can be implemented in a privacy 
preserving fashion, in a manner that accounts for vehicle 
density and composition, and in a manner that is provably 
convergent. Simulations are given to illustrate the efficacy and 



(a) Evolution of the vehicles’ speeds. 



(b) Evolution of the overall energy loss. 

Eig. 11: Simulation results for the network of EVs: Algorithm 
1 is applied until time 0.33 h, and then two different speeds 
(below and above the optimal one) are suggested in [0.33,0.66] 
h and [0.66,1] h, respectively. 


the acceptability of the algorithm. Einally, the algorithm has 
been implemented in a real production vehicle embedded into 
a HIE emulation using nothing more than a smartphone and 
a commercially available OBD-II plug-in. In future work we 
intend to extend the HIE simulation to include more vehicles 
(and thus, more drivers) to better test the compliance of people 
to follow the recommended speed. 


Appendix 

Outline of proof of Theorem[T] 


In this section, we give an outline of the proof for the 
claims in Subsection 


III-A in which we largely rely on the 
results obtained in my 


In Theorem [T] statement (i) is a consequence of the Banach 
contraction theorem. It is a straightforward calculation to 
show that the bounds (|8]l ensure that the function h defining 
the Lure system (|^ is in fact a global strict contraction on K. 
Statement (ii) then follows directly from the definition of h: 
if h{y*) = y* then G{y*e) = 0 and by (|^ this is equivalent 
to the optimality condition Q. The global optimality of y*e 
for the optimisation problem (|^ of this optimal point then 
follows as Q is the standard first order necessary condition 
for optimality, and because strict convexity of the cost 
functions implies this condition is also sufficient. Uniqueness 
is a further consequence of strict optimality. 


It therefore remains to show that statement [T] (iii) holds. To 
this end we recall the following two lemmas from p^. 
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Lemma 2 ( p9) ) Let {-P (fc)}j,gfj^ be a sequence of row- 
stochastic matrices. If {y (fc)}^gj^ is a solution of the Lure 
system then {y {k) is a solution of 0 . 

Lemma 3 ( |^)) Let {P be a strongly ergodic se¬ 

quence of row-stochastic matrices, and suppose that G : K" —)■ 
M is continuous and satisfies the following conditions: 

(i) there exists an e > 0 such that G satisfies a Lipschitz 
condition with constant L > 0 on the set 

Bg (E) := {a; G K" : dist {x, E) < e} , 

where dist {x,E) := inf {||a:: — z\\ : z G E} is the dis¬ 
tance of a vector x G M" to the consensus set E 
span{e}; and 

(ii) there exists constants / 3,7 > 0 such that 

\h{y)\ < |y| -7 when \y\ > /3, 

where h{y) = y -\- G (ye). 

Then, every trajectory of 0 is bounded. 

It is easy to see that G as defined in 0 satisifies the 
conditions of Lemma Indeed, G is even globally Lipschitz 
continuous because of the Lipschitz continuity assumption 0 . 
Furthermore, as is a strict contraction on M with fixed point 
y*, we may denote the contraction constant of by 0 < c < 1 
and obtain for any y gM. such that 

\hiy)\ < \hiy)-y*\-G\y*\ < c\y-y*\-G\y*\ 

<c\y\-G ( 1 +c) \y*\ = \y\- ( 1 -c) \y\ + ( 1 +c) \y*\, 

from which it is easy to derive constants (3 and 7 . 

Finally, if every trajectory of 0 is bounded, then every 
trajectory has a nonempty bounded w-limit set. Because of the 
averaging property of stochastic matrices and the assumption 
of uniform strong ergodicity, this w-limit set is a subset of 
the span of e. By part (i) of the theorem, the Lure system 
has a globally asymptotically stable fixed point. Lemma on 
the other hand ensures that on spanje} the trajectories of 0 
and (|^ (multiplied by e) coincide. It follows that restricted to 
spanje}, the optimisation algorithm (|^ has only one w-limit 
set, namely y*e. It then follows from a continuity argument 
that y*e is a globally asymptotically stable fixed point of 0 . 
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